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In this article, we report the influence of the interfacial roughness on the thermal boundary 
conductance between two crystals, using molecular dynamics. We show evidence of a transition 
between two regimes, depending on the interfacial roughness: when the roughness is small, the 
boundary conductance is constant taking values close to the conductance of the corresponding 
planar interface. When the roughness is larger, the conductance becomes larger than the planar 
interface conductance, and the relative increase is found to be close to the increase of the interfacial 
area. The cross-plane conductivity of a superlattice with rough interfaces is found to increase in 
a comparable amount, suggesting that heat transport in superlattices is mainly controlled by the 
boundary conductance. These observations are interpreted using the wave characteristics of the 
energy carriers. We characterize also the effect of the angle of the asperities, and find that the 
boundary conductance displayed by interfaces having steep slopes may become important if the 
lateral period characterizing the interfacial profile is large enough. Finally, we consider the effect of 
the shape of the interfaces, and show that the sinusoidal interface displays the highest conductance, 
because of its large true interfacial area. All these considerations are relevant to the optimization 
of nanoscale interfacial energy transport. 

PACS numbers: 68.35.Ja, 07.05.Tp, 44.10.+i 


I. INTRODUCTION 


The existence of a finite thermal boundary resistance 
between two solids has important practical consequences, 
especially in the transport properties of nanostructured 
materials. When the distance between interfaces be¬ 
comes submicronic, heat transfer is mainly controlled by 
the interfacial phonon transmission, which in turn gov¬ 
erns the thermal boundary resistance. In certain applica¬ 
tions, such as electro-optical modulators^, optical switch¬ 
ing devices^ and pressure sensors^, a low resistance is de¬ 
sired to favor energy flow. In thermoelectric devices on 
the contrary, a large resistance is preferable so as to gen¬ 
erate large barriers for a wide class of phonons^^— . Two 
strategies may be followed in order to tune the value 
of the boundary resistance between two solids. Either, 
the solid/solid interaction is changed through the cou¬ 
pling with a third body, which is usually a self-assembled 
monolayeri^. The other possibility is to modulate the 
interfacial roughness^. This latter direction has been il¬ 
lustrated experimentally through chemical etchingiSrJ^. 
However, a theoretical model describing the effect of the 
interfacial roughness on the thermal boundary conduc¬ 
tance at room temperature is still lackin g^^d^ . Note that 
the role of the interfacial roughness on the Kapitza con¬ 
ductance has been underlined a long time ago, in the 
context of liquid Helium/metal interfaces at very low- 
cryogenic temperature a^^d^ . At these temperatures, the 
phonon coherence length may be comparable with the 
typical heights of the interface, leading to strong phonon 
scattering which is put forward to explain the high val¬ 


ues of the conductance experimentally reported, as com¬ 
pared with the classical acoustic mismatch theory which 
assumes planar interfaces'^. Such considerations have re¬ 
ceived less attention for room temperatures solids, prob¬ 
ably because in this case the phonon coherence length is 
very small. 

Understanding the role of the interfacial roughness has 
also important consequences in the transport properties 
of super lattices, which are good candidates for thermo¬ 
electric conversion materials, thanks to their low thermal 
conductivityi^. Designing superlattices with rough inter¬ 
faces has been recently achieved, opening an avenue for 
reducing the thermal conductivity in the direction per¬ 
pendicular to the interfaces'^. Again, the physical mech¬ 
anisms at play in the heat transport properties of rough 
superlattices have not been elucidated so far. Molecu¬ 
lar dynamics offers a privileged route to understand the 
interaction between the energy carriers in a solid and 
the asperities of the interface^Sri^. In this article, we 
use molecular dynamics to probe interfacial heat transfer 
across model rough interfaces. Because of the difficulty to 
determine the temperature jump across a non-planar in¬ 
terface, we have used transient simulations which enable 
to compute the thermal boundary resistance characteriz¬ 
ing rough interfaces. In section [Hi we describe the struc¬ 
tures used to probe the conductance of rough interfaces. 
In section imi we explain the methodology retained to 
extract the thermal boundary conductance from molec¬ 
ular dynamics. The simulation results are presented and 
discussed in section lYl We first concentrate on model 
interfaces made of isosceles triangles. 

For these model interfaces, we present the results for 
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the thermal conductance as a function of the interfacial 
roughness and interpret the results using a simple acous¬ 
tic model in subsection IIV B1 In the following section, we 
characterize the effect of the angles of the asperities. Fi¬ 
nally, in sub-section IIV D[ we have appraised the effect 
of the interfacial shape. We discuss the consequences of 
this work in the Conclusion. 


II. STRUCTURES AND SAMPLE 
PREPARATION 

We will consider model rough interfaces, constructed 
from two perfect fee Lennard-Jones solids whose interface 
is orientated along the crystallographic [100] direction, as 
represented in figure |T] We introduce some 2D roughness 
in the xz plane, where x and z denote respectively the 
[100] and [001] directions. As we use periodic boundary 
conditions in all spatial directions, the system studied 
is similar to a superlattice. The dimension in the y di¬ 
rection has been fixed to 10 oq, where oq is the lattice 
parameter, while the dimension in the z direction -the 
super lattice period-has been varied between 5 and 40 ag. 
All the atoms of the system interact through a Lennard- 
Jones potential VLj(r') = 4e ((cr/r)^^ — (a/r)^) truncated 
at a distance 2.5tT. A single set of energy e and diame¬ 
ter cr characterizes the interatomic interaction potential. 
As a result, the two solids have the same lattice con¬ 
stant oq. To introduce an acoustic mismatch between the 
two solids, we have considered a difference between the 
masses of the atoms of the two solids, characterized by 
the mass ratio = m 2 /mi. In all the following, we will 
use mr = 2, which has been shown to give an impedance 
ratio typical of the interface between Si and Ge^^. From 
now on, we will use real units where e = 1.67 10“^^ 
J; cr = 3.4 10“^° m and mi = 6.63 10“^® kg, where 
these different values have been chosen to represent solid 
Argon. With this choice of units, the unit of time is 
T = Y^mcr^/e = 2.14 ps, and the unit of interfacial con¬ 
ductance is G = fcB/(TCT^) ~ bQMW/K/m?. The dif¬ 
ferent interfaces have been prepared as follows: first the 
structures have been generated by mapping the space 
with fee structures using the lattice parameter of the fee 
LJ solid at zero temperature^i: ag{T = QK) = 1.5496cr. 
The structures have been equilibrated at the two final fi¬ 
nite temperatures T = 40 K and T = 18 K using a com¬ 
bination of a Berendsen, a Nose Hoover thermostat and 
a barostat at 0 atm2^. The total equilibration time lasts 
one million time steps which corresponds to a total time 
of 4, 28 ns. The equilibrium lattice parameters have been 
found to be: a = 1.579cr at T=40 K and a = 1.5563cr at 
T=18 K. 

In this article, we considered different types of rough 
interfaces, as represented in Fig.[TJ The first one consists 
of triangular shaped interfaces, having a constant angle 
q; = 45 deg with respect to the xy plane and a variable 
height h. In the second type of interface analyzed, we 
keep constant the interfacial height h and we vary the 
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FIG. 1. (Color online) Diagram showing the different pa¬ 
rameters that will be varied, for the triangular shaped in¬ 
terfaces. Schematic representation of the different parame¬ 
ters that have been varied: Top middle (Fig. 3b): interfacial 
height at a constant value of the angle a = 45 deg; Medium 
(fig 3c): angle q at a constant value of the interface height 
h. Bottom (Fig 3d): angle a and interfacial height /i at a 
constant value of the interfacial area A. 


angle a. In the third case, both the angle and the height 
are varied keeping constant the interfacial area A. Fi¬ 
nally the effect of the shape of the interfaces has been 
also appraised considering totally rough interfaces, small 
triangles juxtaposed on triangular interfaces, square and 
wavy shaped interfaces. This analysis covers all the pos¬ 
sible parameters that might be involved in the geometry 
of the interfaces with the aim to quantify their influence 
on phonon interfacial transport. 


III. THERMAL BOUNDARY CONDUCTANCE 
FROM TRANSIENT SIMULATIONS 

In this section, we briefly review the methodology 
adopted to probe the interfacial conductance between 
two solids, using transient out-of-equilibrium simula¬ 
tions. 

Generally speaking, there are different methods to ex¬ 
tract the boundary conductance from molecular sim¬ 
ulations: either a net heat flux g is applied through 
the coupling of two energy reservoirs, and one mea¬ 
sures the finite temperature jump AT across the inter- 
fac o^^i^^ . This allows to compute the interfacial conduc- 
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FIG. 2. (Color online) Energy decay of the heated solid ob¬ 
tained using transient non-equilibrium simulations. Dashed 
lines show the exponential £t. The parameters are: Total 
Length=40 oq; temperature T = 40 K; mass ratio nir = 2. 


identical solids, and calculated the corresponding internal 
resistance. The Kapitza conductance calculated in this 
article has been obtained after having substracted this 
internal resistance: 

1/Gk = 1 /Gi 2 - 1/Gn (3) 

where G12 is the conductance measured for the interface 
between solid 1 and solid 2; and Gn is the conductance 
measured between identical solids using eq. [T] This pro¬ 
cedure has been followed for all the systems studied in 
this article. Finally, for the simulations discussed in this 
article, we have used between 5 and 10 independent con¬ 
figurations depending on the system size, to determine 
the value of the Kapitza conductance, and the error bar 
has been found to be typically 15 percents. 

IV. RESULTS 


tance Gk = 9/AT. For the rough interfaces that we will 
consider in the following, it may be difficult to clearly 
identify a temperature jump, especially if the roughness 
is large. On the other hand, transient non-equilibrium 
simulations do not require to resolve spatially the tem¬ 
perature field in the vicinity of the interface, and for this 
reason are well adapted to the determination of the con¬ 
ductance of imperfect interfaces. The principle is akin to 
the thermoreflectance technique, and consists in heating 
instantaneously one of the two solids and record the tem¬ 
poral evolution of the temperature of the hot soli d^^i^^ . 
The conductance Gk is then obtained from the time r 
characterizing the thermal relaxation of the hot solid: 


_ miks 

AAt 


( 1 ) 


where Ni is the number of atoms of solid 1, ks is the 
Boltzmann’s constant, and A is the interfacial area. Be¬ 
cause the temperature of the heated solid may display 
some oscillations which may make the determination of 
the time constant r difficult, we have rather used the 
decay of the energy Ei to extract r: 

iGi j.fcei 


where the second term corresponds to the interatomic 
potential, here supposed to be pair-wise. An example of 
the time decay of the energy during thermal relaxation 
of the hot solid is displayed in Fig. [21 showing that the 
exponential decay hypothesis is reasonable. 

In practice after equilibration of the system, we have 
heated one of the two solids by an amount of 18 K and 
followed the thermal relaxation of the hot solid during 
a time interval after which the energy has decreased by 
a factor 10. To remove any contribution stemming from 
internal phonon scattering in the heated solid, we have 
run in parallel simulations across the interface between 


In this section, we present the simulation results ob¬ 
tained using the transient simulations, as detailed above. 
We will successively study the effect of the interfacial 
roughness, the angle of the asperities, and the shape of 
the interface. A summary of the different parameters 
that will be varied is depicted in Fig. |T] and Fig. |9l 


A. Effect of the superlattice period and number of 
periods 

In this section, we first quantify finite size effects in 
the determination of the conductance of rough interfaces, 
as measured by eq. [1] It is important to note that the 
system simulated is not a single isolated interface, but 
rather a superlattice because of the periodic boundary 
conditions. It is thus relevant to assess the influence of 
the superlattice period on the thermal conductance as 
measured by eq. [H To this end, we will consider model 
rough interfaces, made of isosceles triangles, as depicted 
in fig.dl In figure m we report the conductance of trian¬ 
gular shaped interfaces having a fixed roughness height 
32 MLs, and a varying period. For the two temperatures 
considered, the thermal conductance is found to decrease 
with the system size for short periods, and then satu¬ 
rates for periods larger than 30 nm. The increase of the 
conductance for thin layers may be explained because 
long wavelength phonons will not see two independent 
interfaces but rather a single one. A similar trend has 
been also reported in lattice dynamics simulations^ and 
Green function calculations^!. For thick layers, the con¬ 
ductance measured is constant, and has converged to the 
value characterizing an infinitely thick film. 

We remark also that the conductance is higher at high 
temperatures. Generally speaking, the thermal bound¬ 
ary conductance is found to increase with temperature, a 
trend often attributed to the existence of inelastic phonon 
scattering at the interface!^!!!. This behavior is consis¬ 
tent with our simulation data. In the following, we will 
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FIG. 3. (Color online) Thermal boundary conductance of 
isosceles triangular shaped interfaces having a roughness 
height h — 32 MLs, as a function of the period p defined 
in fig[T] Top: T = 40 K; Bottom: T = 18 K. Lines are guides 
to the eye. The mass ratio is rtir = 2. 

fix the period to 20 ag as it leads to moderate finite size 
effects as already found for superlattices with planar in¬ 
terfaces^. Finally, since the system we simulate is akin 
to a superlattice because of the periodic boundary con¬ 
ditions, it is important to probe the effect of the number 
of periods on the measured conductance. Figure 2] quan¬ 
tifies the effect of the number of interfaces on the con¬ 
ductance measured in transient simulations. From this 
figure, we can conclude that within error bars, the num¬ 
ber of interfaces has a mild effect on the conductance 
that we calculate. This is the behaviour expected, as 
we probe a quantity characterizing the interface solely, 
independently on the number of interfaces. 

B. Effect of the interfacial roughness 

We will now concentrate on the influence of the inter¬ 
facial roughness on the thermal boundary conductance. 
We will consider rough interfaces having an angle a fixed 
at a = 45 deg, while the height h of the interface is in¬ 
creased so as to change the interfacial roughness, as seen 
in fig [T] b. For the following, it is important to keep in 
mind that, when varying the height of the interface h at 
a constant value of the angle, the true interfacial area 
remains constant, and larger by a factor 1/ cos a = y/2 
than its corresponding projection on the horizontal xy 
plane. 

Figure 5 displays the evolution of the measured 
Kapitza conductance, as a function of the interfacial 
roughness for two temperatures. The conductance of a 
planar interface, which corresponds to the value h = 0 
has been also indicated for the sake of comparison. Two 



FIG. 4. (Color online) Thermal boundary conductance of 
isosceles triangular shaped interfaces having a roughness 
height h = 32 MLs and a Period 40 ag, as a function of 
the number of periods of the superlattice. Lines are guides to 
the eye. The temperature is T = 40 K and the mass ratio is 
nir ~ 2. 


regimes are to be distinguished, depending on the rough¬ 
ness of the interface h. When the height is smaller than 
typically 20 monolayers, the conductance seems to be 
constant, or slightly decreases with the roughness, taking 
values close to the planar interface conductance. When 
the interfacial height becomes larger, the conductance 
suddenly increases and tends to saturate for very rough 
surfaces. 

Interestingly, the increase of the conductance between 
planar and very rough surfaces is found to be close to 
the increase of the interfacial area. This is materialized 
in fig. 5, where we have shown with dashed lines the 
value of the conductance obtained by multiplying the 
conductance of a planar interface by a factor 1/cosa. 
Depending on the superlattice period, the increase of the 
Cross-Plane thermal conductivity has been found to be 
included between 1.3 and 1.5, which encompasses the 
value ~ 1.41. This reinforces the message accord¬ 
ing to which the thermal conductivity of a superlattice 
is mainly controlled by the Kapitza resistance exhibited 
by the interfaces, which in turn seems to be primarily 
governed by the interfacial area. 

We now give some qualitative elements to understand 
the previous simulation results, regarding the influence 
of the interfacial roughness on the Kapitza conductance. 
At this point, it is important to have in mind that in 
the situations that we have modeled, the energy carri¬ 
ers are phonons which are classically populated. A given 
phonon mode is characterized among others by its wave¬ 
length A, which may take practically any value between 
an atomic distance 2a and the simulation box length L^. 
First, let us concentrate on the case of a small roughness 
h, as exemplified in fig. [51 In this case, the majority 
of phonon modes have a wavelength larger than h, and 
they see the interface as a planar one: the transmitted 
heat flux is then controlled by the projected area. On 
the other hand, when the interface is very rough, most of 
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FIG. 5. (Color online) Thermal boundary conductance of 
isosceles triangular shaped interfaces, as a function of the in¬ 
terfacial roughness. Top: T = 40 K; Bottom: T = 18 K. 
We have also indicated the conductance of the corresponding 
planar interfaces (ML=0) and very rough interface (ML=60). 
The horizontal dashed lines show the conductance obtained 
after rescaling the conductance of the planar interface by the 
true interfacial area. The solid lines denote the theoretical 
model eq. m, with the parameter ^ = 0.2. The parameters 
are: total Length=40 oq; mass ratio = 2. 


the phonons have a wavelength smaller than the height 
h, obviously the phonons no longer feel the interface as 
planar, phonon scattering becomes completely incoher¬ 
ent, and the transmitted heat flux is controlled by the 
true surface area. 

To put these arguments on quantitative grounds, we 
will consider the following expression of the thermal con¬ 
ductance, inspired by the classical AMM model^^. We in¬ 
troduce a mode-dependent fraction ipWi which depends 
on the considered wavelength, and which is equal to 1 
when the wavelength is supposed to be small compared 
with the interfacial roughness h, and equal 0 in the op¬ 
posite case. We define a dimensionless parameter such 
that: 


tp{X) = 1 if A < 

ip{X) = 0 otherwise (4) 


The parameter ^ will be the adjustable fitting param¬ 
eter of the model. The interfacial conductance is then 
supposed to be given by : 

Gk = T:CpkBCi{ / g{Lo)il){X)dujIi2 
^ Jo 

+ ^/ g(w)(l -'0(A))du;Ji2) (5) 

where p is the crystal number density, ci is the average 
sound velocity in medium 1, womin is the Debye frequency 
of the softer solid. The parameter C is a scaling factor 
which accounts for the tendency of the AMM model to 
overpredict the measured Kapitza conductance^^. The 
integral I12 involves the angular dependent transmission 
coefficient: 


Ti2 = 


'iZiPiZ2P2 , 

0 {Zipi + Z2P2)^^^ 


( 6 ) 


where Zi = p^Ci are the acoustic impedances of the two 
solids, p™ being the mass density, and pi = cos^i is a 
shorthand notation to denote the cosine of the phonon 
incident angle^^. Finally, the quantity AjAo is the ratio 
of the true interfacial area over the projected one. The 
physical motivation of eq. is simple: phonons hav¬ 
ing a wavelength A larger than contribute to a trans¬ 
mitted heat flux proportional to the projected area Mo, 
while phonons modes having a wavelength smaller than 
contribute to the transmitted heat flux proportionally 
to the true surface area. We have compared the predic¬ 
tion of eq. [5] to the simulation results discussed before. 
To this end, we have assumed Debye solids, with a vibra¬ 
tional density of states g{(jj) = a;^/(27r^cf), and for the 
sake of consistency, the mode-dependent wavelength A 
has been taken to be simply related to the frequency w : 
A = 2ttC l joj. Figure 5 compares the predictions of eq. [5] 
to the simulations results. The values of the planar inter¬ 
face conductance has been rescaled by a factor C = 3 and 
4 at the temperatures T = 18 K and T = 40 K respec¬ 
tively. These correction factors account for the fact that 
the simple AMM model relies on several assumptions- 
Debye solids-elastic scattering-which may lead to a dis¬ 
crepancy with the MD value. We have chosen the two 
factors because they yield good agreement with the MD 
value for smooth interfaces. Apart from this rescaling, 
the parameter ^ has been treated as the only fitting pa¬ 
rameter. Figure 5 shows that a good agreement is found 
using the value ^ = 0.2, for the two temperatures con¬ 
sidered. The smallness of the fitting coefficient may be 
understood in the following way: consider a given phonon 
mode: if its wavelength is larger than the roughness h, 
the effective scattering area would be the projected one. 
On increasing the roughness, h will become comparable 
with A, and the interface will strongly scatter the con¬ 
sidered phonon in all directions. This will contribute 
to slightly diminish the conductance, as compared with 
the planar case, in agreement with the simulation data 
points. It is only when the roughness becomes very large 
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FIG. 6. (Color online) Schematic representation of the rough¬ 
ness indnced phonon scattering. Top: case of a small rough¬ 
ness. The huge majority of incoming phonons see the interface 
as a plane, and the transmitted heat flux is proportional to the 
projected interface area. Bottom: case of a large roughness. 
Most of the phonons have a wavelength larger than the inter¬ 
facial roughness, and the transmitted heat flux is proportional 
to the true surface area. For the sake of the representation, 
we have not drawn the reflected waves. Note also that the 
phonon wavelength is generally not conserved at the passage 
of the interface. 

as compared with the wavelength A <C h, that interfacial 
scattering becomes again negligible and the transmitted 
energy is proportional to the true area. The fitting pro¬ 
cedure concludes that this regime is reached when the 
wavelength becoirres smaller thair typically one fifth of 
the interfacial roughiress. 


C. Effect of the angle of the asperities 

So far, we have considered the case where the angle 
a was constant. We now discuss the effect of varying 
the slope of the model interfaces on the interfacial en¬ 
ergy transfer. First, we will change the angle at a fixed 
value of the interfacial height h, as represented in fig. [T] 
c. Figure o shows the evolution of the Kapitza conduc¬ 
tance as a function of the angle, at the two temperatures 
considered. The constant height h retained here corre¬ 
sponds to the regime of large roughnesses in terms of 
the figure 5 discussed before. We have also indicated the 
conductance of a planar interface, for the sake of compar¬ 
ison. The evolution of the conductance with the asperi¬ 
ties angle seems to be non-monotonous: first, it increases 
for low angles, reaches a maximum for an asperities an- 



FIG. 7. (Golor online) Thermal boundary conductance as 
a function of the angle of the asperities a. The height of 
the asperities is fixed here and equal to 24 ML. The points 
corresponding to a = 0 denote the conductance of a planar 
interface. The solid lines show the interfacial conductance 
rescaled by the true surface area: Gk = Gk^ch = Q)/ cos a. 
The other parameters are: Total length=40 no; mass ratio 
rrir = 2. 



angle (degrees) 


FIG. 8. (Golor online) Same as fig. [T] but for a constant 
value of the true surface area. The interfacial heights are re¬ 
spectively h = 21,34,43 and 45 monolayers for the asperities 
angles a = 25.8,45,64 and 71.7 deg. The solid lines show 
the interfacial conductance rescaled by the true surface area: 
Gk = Gk {Q- = 0)/ cos a. 


gle between 30 and 45 degrees, and then decreases when 
the angle becomes large. In particular, the conductance 
for an angle greater than 60 deg becomes smaller than 
the planar interface conductance. This is all the more 
remarkable as in this latter case, the true surface area 
may increase by a factor 4 as compared to the planar 
interface. This discrepancy is best shown after compar¬ 
ing the simulation results with the rescaled conductance 
Gif (a = 0 )/cosq;, which accounts for the increased sur¬ 
face area induced by the asperities. It is immediately 
clear, that for the lowest values of the asperities an¬ 
gles, a = 30 and 45 deg, the rescaled conductance seems 
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random rough squares wavy 


FIG. 9. (Color online) Illustration of the different interfacial 
shapes simulated, namely random, rough, square and wave¬ 
like interfaces. 


to describe reasonably interfacial energy transfer. On 
the other hand, at large values of a, the theoretical ex¬ 
pression overestimates greatly interfacial transport. Two 
phenomena may explain the poor conductance reported: 
first, on increasing the angle, phonon multiple scattering 
and back-scattering may contribute to diminish interfa¬ 
cial transmission. This has been evidenced by Rajabpour 
et al. using Monte-Carlo ray tracing calculations^!. Sec¬ 
ondly, for the steep slopes interfaces considered here, the 
effective surface area seems to be the projected one, not 
the true area, even if the height of the asperities is large. 
This may be understood qualitatively because for steep 
interfaces, even if the height is large, the lateral corre¬ 
lation length I = h/tana may become comparable with 
the phonon wavelengths, and the effective interfacial area 
becomes the projected one. For these steep interfaces, the 
regime where the transmitted heat flux is controlled by 
the true surface area should occur at a very large value 
of the interfacial height h. To verify this assessment, we 
have run simulations where the true surface area has been 
kept constant (cf fig. [T]d). The results are displayed in 
fig. [51 which concludes a different scenario as compared 
to the evolution shown previously in fig. 0. The evo¬ 
lution of the conductance with the angle is no longer 
non-monotonous as previously observed, but it rather in¬ 
creases monotonously with a. For the relatively small 
values of the angles a, the conductance measured may 
even exceed the rescaled one. We have no interpretation 
for these large values reported here. Increasing further 
the angle a, the simulation data takes values close to the 
scaled conductances Gk{<^ = Oj/cosa. Note in particu¬ 
lar, that the increase of the conductance is pretty large, 
overpassing the conductance of a planar interface by a 
factor larger than 3. In this regime, and for these steep 
interfaces, it is highly probable that the regime of rough 
interfaces, in the terms of the previous discussion has 
been reached: heat transmission becomes controlled by 
the true surface area. These large enhancement of the 
Kapitza conductance open the way to design interfaces 
with tailored interfacial energy transport properties. 



FIG. 10. (Color online) Thermal boundary conductance for 
the different interfacial shapes represented in fig. (9) The 
height of the different interfaces is fixed here and equal to 
12 ML. The other parameters are: Total Length=40 uq; mass 
ratio nir = 2. 


D. Effect of the shape of the interfaces 

We end up the presentation of the results by discussing 
the effect of the shape of the interface on the boundary 
conductance. All the previous discussion concentrated on 
model triangular interfaces, and it is worth asking how 
general are the conclusions drawn from the study of the 
particular type of surfaces. To appraise this question, 
we have considered different shapes of the interfaces, as 
depicted in fig. (9] The common characteristics of these 
surfaces is the mean interfacial height, here fixed at a 
value h = 12 MLs. Different morphologies have been 
designed, ranging from the random interface, to the case 
of square-like surface and wavy interface obtained by a 
sinusoidal modulation of the interfacial height. 

Figure do] compares the interfacial conductance for the 
different shapes shown before, at two different tempera¬ 
tures. The relatively large values reported at the highest 
temperature may be explained by inelastic phonon scat¬ 
tering taking place between two interfaces. The shape 
of the interface seems to affect considerably interfacial 
transport: random interfaces display a conductance prac¬ 
tically equal to the planar interface. Rough interfaces 
may transfer energy slightly better than planar interfaces 
depending on the temperature. On the other hand, wavy 
and square-like interfaces tend to favor energy transmis¬ 
sion, the wavy pattern displaying the highest conduc¬ 
tance among the different shapes analyzed. These re¬ 
sults may be interpreted qualitatively: random and rough 
shaped interfaces display a distribution of length scales, 
which tend to promote phonon scattering: even if the 
global height h is large, in terms of triangular shaped 
interface, the effective area for the phonons is not the 
true surface area, but rather the projected area, because 
h is not the only relevant roughness parameter, and the 
interaction between incident phonons and small length 
scale asperities tend to diminish the effective transmis- 
















sion area. On the other hand, regular shaped pattern 
do not display such a distribution of length scales, and 
interfacial heat transport becomes controlled by the true 
surface area: as soon as the majority of phonon modes 
has a wavelength greater than the single length h charac¬ 
terizing the interfacial morphology, one enters in a ” large 
roughness” regime where the energy transport becomes 
governed by the true surface area, and the conductance 
is increased as compared with the planar case. The in¬ 
terfacial conductance is found to be the highest for the 
wavy interface, because it has the greatest surface area. 

V. CONCLUSION 

In summary, we have concentrated on the role of the in¬ 
terfacial roughness on energy transmission between solid 
dielectrics. Thanks to the versatility of molecular dy¬ 
namics simulations technique, we have conceived model 
rough surfaces, and probe their ability to conduct heat. 
The scenario emerging from the simulations is the follow¬ 
ing: when the roughness introduced is small, most of the 
phonons see the interface as a planar one and the effective 
surface area contributing to the transmitted heat flux is 
the projected area, not the true one. In this regime, one 
does not expect a Kapitza conductance much different 
from the planar interface. On the other hand, when the 
roughness becomes large enough, typically 20 monolay¬ 
ers in our case, most of the phonons propagating towards 
the interface are incoherently scattered, and the effective 
surface area becomes the true surface area. This latter 
may differ significantly from the projected one, and this 
is the reason why the boundary conductance of rough 
interfaces may be greatly enhanced, as compared to pla¬ 
nar interfaces. This has been demonstrated in this work, 
with the example of triangular shaped interfaces display¬ 
ing steep slopes: provided the lateral dimensions char¬ 
acterizing the interfacial roughness is large enough, the 
increase of the conductance may be threefold. On the 
other hand, we have probed the conductance of randomly 
rough interfaces, and shown that they display in general 
conductances comparable or smaller than atomic planar 
interfaces. This difference of behavior is explained by the 
distribution of length scales displayed by the randomly 
rough surfaces, in comparison with our model patterned 
surfaces. 

The roughness analyzed in this article was large com¬ 
pared to the lattice constants. The case of atomic rough¬ 
ness has been more widely addressed in the literature, 
and wave-packet simulations^^ give a clear picture of the 
effect of small atomic roughness on phonon transmission: 
long wavelength phonons see the interface as ideal and 
do not contribute to the change of the thermal bound¬ 
ary conductance. On the other hand, short wavelength 
phonons strongly interact with the small scale roughness. 


and the corresponding change in phonon transmission is 
found to depend on the structure of the interface: for reg¬ 
ular shaped patterned interfaces, constructive wave inter¬ 
ference lead to enhanced transmission thereby increasing 
the boundary conductance^. Random atomic rough¬ 
ness promotes incoherent phonon scattering, reducing the 
thermal conductance. These observations are consistent 
with MD results concerning the cross-plane conductivity 
of super lattices with rough interfaces: for regularly pat¬ 
terned interfaces, the cross-plane conductivity is slightly 
greater than ideal superlattices^S, and the boundary con¬ 
ductance is enhanced^. When the roughness is random, 
the cross-plane conductivity shows a small decrease as 
compared with planar interface a^^i^^ . The small ampli¬ 
tude of the reduction is related to the small proportion of 
energy carriers affected by the atomic roughness. Small 
reductions have been also reported for Si/Ge superlat¬ 
tices with a one layer of interfacial mixing in the incoher¬ 
ent regime of transport^. If it is reasonable to rationalize 
such variations in terms of an atomic interfacial rough¬ 
ness, it is less clear for super lattices having thicker mixed 
layer. Large enhancements have been observed in this 
latter case, using MD^S. Further work is clearly needed 
to understand if part of these enhancements is explained 
by the large scale interfacial roughness^i 

Most of the results reported here concern regular 
shaped patterned interfaces. MD results seem to con¬ 
clude that these patterned interfaces are good candi¬ 
dates to enhance the intrinsic boundary conductance be¬ 
tween two semi-conductors. On the other hand, ran¬ 
domly rough surfaces should be considered if one prefers 
to reduce the Kapitza conductance between two solidsi^. 
In particular, in the context of superlattices, randomly 
rough interfaces should be designed if one aims at tai¬ 
loring materials with the lowest cross-plane thermal con¬ 
ductivity. 

We have also introduced a simple model to rationalize 
the variations of the thermal boundary conductance as 
a function of the interfacial height of our model rough 
interfaces. Further analytical work is clearly desired to 
understand the interplay between the interface morphol¬ 
ogy and energy interfacial transport. This will enable to 
define new directions for the design of interfaces with op¬ 
timized energy transport properties, with a relative low 
cost. 
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